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We present theoretical investigation of a single magnetic impurity in a d-wave superconductor 
using the large N limit. It is shown that the Kondo screening occurs only in the presence of 
the particle-hole asymmetry. We find analytical expressions for the Kondo temperature, magnetic 
susceptibility and scattering matrix near the phase transition. The results are generalized for the 
density of states vanishing with an arbitrary exponent. Also we briefly study the modifications of 
the theory for the case of a non-magnetic impurity which induces a staggered spin configuration on 
the nearby copper atoms. 



h 

CO 



I 

O 

o 



> 

in 

oo 

o 



T3 

o 



X 



I. INTRODUCTION 

Recent experiments in d-wave BSCCO superconduc- 
tors showed some interesting phenomena in the presence 
of impurities on copper sites. In particular, STM showed 
a strong variation of the electron density of states (DOS) 
near Ni jj], |2) and Zn atoms in E^S^CaC^Os+a- Us- 
ing NMR Kondo- like effect was in cuprates doped by Zn, 
Ni or Li impurities both above and below T c [§, [|, |[ @, f| , 
and an uncompensated spin S — 1/2 appeared near the 
Zn impurity ||. Although it is tempting to explain 
the tunneling experiments by the pure potential scatter- 
ing (n]] , there are several key features which are not easily 
justified within this model (see e.g. jll)). For example, 
the low frequency resonance can be achieved only for a 
large non- universal value of the potential coupling p2| . 
On the other hand, the Kondo screening gives a natural 
low energy scale and the results appear to be much more 
robust with respect to the variation of the composition 
and coupling constants |ll|. We note that using STM 
data alone it is not possible to distinguish between the 
potential and Kondo scattering mechanisms leading to 
the variations of the DOS . The screened impurity acts 
as a spin singlet and is therefore similar to the poten- 
tial scatterrer. However the unitarity limit in the Kondo 
case occurs at the coupling slightly exceeding the critical 
value, which is not necessarily large. The other important 
advantage of the Kondo mechanism is that it simultane- 
ously explains NMR and STM measurements. 

The cuprates are unconventional superconductors with 
linearly vanishing density of states near the four nodes. 
The usual Kondo picture is not valid there jlj, [lj, [H], 
[I?], [l^, |l9|]. In particular, there is no phase transition 
if the spin-spin coupling is relatively small so the impu- 
rity remains effectively decoupled from quasiparticles. If 
the DOS vanishes at the Fermi surface slower than the 
square root of energy: p(e) oc |e| r with r < 1/2, then it 
is possible to develop a perturbation theory UM, which 



predicts Kondo transition at the spin-spin coupling larger 
than a certain critical value proportional to r. However 
if r > 1/2 this simple picture doesn't work. Thus both 
numerical renormalization group (NRG) |2(|, pl"[ ] and dy- 
namic large N multichannel approach |22)| showed that 
the Kondo phase exists only if the Hamiltonian doesn't 
respect the particle-hole symmetry. In the symmetric 
case the spin-spin coupling always flows to zero. On the 
other hand, in |2J] it was argued that if the particle-hole 
symmetry is preserved, then there is no screening at the 
level of the total impurity susceptibility T\imp — * 1/4 at 
T — > 0, but still there is a quenching of the impurity spin 
and Txioc - * in the zero temperature limit. 

The model Hamiltonian for a pure superconductor to 
be exploited in the present paper is: 



H = (e k r« + A k r«) ^(k), 



(1) 



where £t = — £o(cos k x +cos k y ) — fi and A k = A (cosA: x — 
cosfc y ) are the kinetic energy and superconducting gap 

respectively, ip(k)=(c\ (k), c{(— k)) is the Gor'kov-Nambu 
spinor. The interaction with a single impurity can be 
expressed in the Coqblin-Schrieffer form |2j| : 

H int = Ctyt( ro )T«Vi(ro) 

-^{(r )F T - FWr ))(F T V T (r ) - ^(ro)^), (2) 



where U and J are the potential and spin-spin cou- 
pling constants, F is the effective impurity spinor in the 
Nambu representation (it is related to the creation and 
annihilation operators by Ff = /f, -Fj. = /f). The single 
impurity occupancy corresponds to the condition: 



flfi = 1 



?t_ij p. _ 



Fj = 0. 



(3) 
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In fact, there is an additional potential like scattering 
term proportional to J in (|^), however, it disappears 
in the path-integral formulation of the problem [^6|. If 
fj, = then the Hamiltonian (|l|) respects the particle 
hole symmetry, i.e. under the transformation k x — > 
k x + 7T, k y — > k y + 7T both the kinetic energy Sk and 
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the order parameter change their signs. This addi- 
tional symmetry is not generic, it disappears if cither 
chemical potential or second neighbor hopping becomes 
not zero. However this notion is very useful, since in 
many cases the particle-hole asymmetry can be treated 
as a perturbation. The model above but without the 
potential scattering term was already examined by C. 
Cassanello and E. Fradkin in [l7). The authors con- 
cluded that there is a Kondo phase transition above some 
critical spin-spin coupling J c . They found that for the 
linear DOS, the Kondo temperature vanishes exponen- 
tially near the phase transition. Also it was argued that 
in the Kondo phase the impurity susceptibility vanishes 
at T — > as Tlnl/T. However, these results were ob- 
tained even in the particle-hole symmetric case. We find 
similar dependences of the Kondo temperature and sus- 
ceptibility only if the particle-hole symmetry is broken, 
e.g. by non-zero potential scattering. The reason in such 
a discrepancy is that the symmetrization of the action 
with respect to positive and negative frequencies used 
in (if], [Tt| is not valid, which resulted in essential incon- 
sistencies there. We come back to this point in the next 
section. The other difference between our approaches is 
that we start from the Hamiltonian ([!]) instead of its lin- 
earized nodal version (if], |l7|. This is important for the 
analysis because (i) the particle-hole asymmetry (U 7^ 
or fj, 7^ 0) is essentially a high energy effect, (ii) the spa- 
tial structure of the scattering can not be reproduced 
from the nodal version of the Hamiltonian, (iii) the ex- 
pressions for the T-matrix and Green functions are much 
more apparent in the original Gor'kov-Nambu represen- 
tation than in that of the nodal particles. Although at 
the later stage we use linearized version of the Hamilto- 
nian for the sake of simplicity and analyticity, the results 
can be easily generalized for the arbitrary spectrum of 
quasiparticles Q. 

We also study in detail the phase diagram at finite 
temperature, where the effects of particle hole asymme- 
try become crucial. The model is generalized to the case 
of an arbitrary exponent r larger than 1/2, where the 
simple Withoff-Fradkin picture |L4| is invalid. We would 
also mention that the problem addressed here but for the 
Anderson magnetic impurity was studied in ||l8| using the 
slave-boson mean field approach. The impurity potential 
had been chosen in such a way that only zero or single 
occupancy was allowed. As a result the Hamiltonian was 
strongly particle-hole asymmetric from the beginning. It 
was found that this model predicts the Kondo effect and 
the particle hole symmetry of the pure superconductor is 
not important. The model also results in a low-energy 
resonance in the DOS and is consistent with the experi- 
mental spatial dependence. 

The paper is organized as follows. In Section II we 
solve the saddle point equations and study the phase di- 
agram in the space of the potential and spin-spin cou- 
plings U and J. It is found that the transition to the 
Kondo phase occurs only in the particle-hole asymmetric 
case. At zero temperature the critical spin-spin coupling 



J c saturates as U — > if r < 1 and diverges as l/y/U if 
r > 1. In the case of a d-wave superconductor (r = 1), 
the Kondo temperature is found to vanish exponentially 
near J = J c in agreement with fl^j . At the fixed temper- 
ature and for r < 1 the crossover from the high temper- 
ature (T > T K ) to the low temperature (T < T K ) limit 
occurs at J diverging as U — ► 0. 

In Section III we examine the expressions for the local 
and total magnetic susceptibilities of the impurity and 
the quasiparticle density of states. It is shown that in the 
particle-hole symmetric case for r < 1 and J > J c the 
local susceptibility diverges as T — » 0, but T\ioc — > 0, 
while T\imp — * const. In the Kondo phase xioc satu- 
rates and Ximp — > 0. Both xioc and Ximp as functions of 
temperature have a maximum at T « Tk- In the end 
of the section we derive the expression for the induced 
quasiparticle DOS near the impurity. It is shown that 
in the Kondo phase the effective scattering matrix has 
a pole at the frequency close to the Kondo temperature 
(f2 w Tk). The crucial difference between this model and 
simple potential scattering | flo| is that the unitarity limit 
here occurs at finite magnitude of the spin-spin coupling 
J ps J c , which doesn't have to be large. 

Section IV is devoted to the discussion of the modifi- 
cations of the model for the case of a non-magnetic im- 
purity, which induces a staggered magnetization nearby. 
This study is motivated by the NMR experiments j|, ||] 
unambiguously showing existence of the magnetic mo- 
ments in the vicinity of the Zn atom, which substitutes 
Cu. From the theoretical point of view this issue is inter- 
esting, since the spin-spin interaction becomes non-local. 
Assuming that the effective spin is induced on the first 
Cu neighbors of Zn we show that the action has four dif- 
ferent saddle points, and for the whole range of coupling 
constants the configuration with the D symmetry of the 
order parameter corresponds to the lowest energy. The 
critical spin-spin coupling for the non-local spin is found 
to decrease with U as opposed to the local spin case. The 
spatial dependence of the induced DOS is also different. 
It is characteristic for the scattering on the four rather 
than the single site and gives an excellent agreement with 
the STM experiments p. 



II. SADDLE POINT APPROXIMATION FOR A 
SINGLE MAGNETIC IMPURITY 

After the standard Hubbard Stratonovich decoupling 
and integrating out Fermion fields the partition function 
for the system with the Hamiltonian given by (|l|) and (|^) 
becomes a functional integral over the auxiliary fields <fi 
and e with the action: 

S= - Ti]n(J- + e{r)r t -<l,(T)G(jAif>(r)\ 


+ ^{rfdr, (4) 
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where trace is taken over all antiperiodic functions of r 
(/(r + /3) = -/(t)); G(5/ar) = G(r , r , 0/flr) is the 



J 



quasiparticle Green function in the presence of the po- 
tential scattering on the impurity: 



G(r, r', u n ) = G°(r - r', w„) - G°(r - r , w B )tfr, (l + Ut z G°{0, w„)) G°(r - r', w„), 

I 



(5) 



and G° is the free electron Green function, which in the 
momentum space is equal to: 



GkM = (- iuJ + ^kT z + A k T x 



(6) 



The Lagrange multiplier e in (|f|) enforces the single oc- 
cupancy constraint. In general is a complex field, how- 
ever its phase can be reabsorbed into e by a simple gauge 
transformation [ g3j . 

It is not hard to show that in a d-wave superconductor 
the free particle Green function at r = is expressed 
as C0|: 



G°(0,w) = G o o (lj) + G o 1 (lj)t z 



(7) 



with Gq and Gj being even and odd functions of fre- 
quency, respectively. If the free Hamiltonian is also 
invariant under the particle-hole transformation then 
G\(lu) = 0. For simplicity we assume this is the case. 
In fact, it can be shown that the models with U ^ and 
Gl can be mapped to each other. Explicitly Gq(oj„) 
is given by 



G8K) 



E 

k 



(8) 



At small frequencies the main contribution to G[j comes 
from the wavevectors close to the nodes, so we can do 
the summation over k and obtain: 



G°K)«^ln 



iu n , loI + A 2 



^7, 



(9) 



where v = y/vpVA is the mean geometrical velocity near 
the nodes and A denoting the upper cutoff of the order 
of A. For the Hamiltonian ([!]) with /i = 0, vf = Eq\/2, 

V A = A V2. 

From the equation (0) we see that the action splits 
into two identical terms for the spin up and spin down 
polarizations so that: 

S= iVTrln(!- + e -0t (r)GT (J_^ 0(r) 

P 

+ j J dT\<t>{r)\*-^ef3, 


where N — 2 is the number of the spin channels, 
G T (l\=G% (I) (l + UG^ 9 "^ 1 



dr 



\dr 



Or 



(10) 



(11) 



Note that there is an additional term "—N/2 e/3" in (|1C 
which enforces the symmetry e — ► — e in the absence of 
the potential scattering. In (10) we find it more conve- 
nient to use a gauge, where <j) is complex, but e is time 
independent. In the large N limit the action is given by 
the saddle point approximation J23|. So it is necessary 
to find the stationary point with respect to the auxiliary 
fields e and d>. Let us start with e: 



2 dS_ 



-G T K)|# 



(12) 



where S is an infinitesimal parameter showing the correct 
procedure of closing the Wick's contour. If U = then 
for any </>: 



I(e)=(3(-l + 2g F (e)) 



J(e)+I(-e) = 0, 



(13) 



where <?f is the Fermi function. In fact this identity 
is just a consequence of the particle-hole symmetry. If 
U 7^ the relation above generalizes to 



I(e,U)+I(-e,-U)=0. 



(14) 



So in the absence of the potential scattering, e = at 
the saddle point. With help of (13), we write the saddle 
point condition as follows: 



E 
E 



(iw„-e+H 2 G T (w„)) (iuj n + \<j>\ 2 G° (uj n )) 
f/|0| 2 G°(a;„)G T (a;„) 



(iLUr, 



-|</>| 2 G T (w„)) (itOr, 



2 G[]K))- 



(15) 



so that both sums become convergent at w„ — > 00. 
Clearly e = as long as \cf>\ = 0. Near the phase tran- 
sition the order parameter \<f)\ is small and therefore e is 
also small. After introducing dimensionless parameters 



/3->A/3, U- 



AU , AJ 2 |0| 2 



J- 



A ' ' rrr 

and using £< 1 we can simplify (p5|) to 



(16) 



tanh ■ 



e/3 



1 



2H 2 ln^ 



2 + 4H 2 lni ) F(o.f i. (17) 



where 
F(0,C7) 



f/.i: 



C/|0| 2 



[ln-^l + ^ + l^ 



f/ 2 a; 2 



(18) 



4 



The function F increases linearly with U at small U and 
saturates at 1/2 for large U. At finite temperature and 
small U or \<f)\, ( |l7| ) shows that e oc T up to logarithmic 
corrections, i.e. instead of ( |l7| ) we can write: 



tanh 



e(3 



l + 2|(/.| 2 ln/3 



2(l + 2|0| 2 ln/3) J F(0,C/). (19) 



As long as T is not too small we can ignore weak loga- 
rithmic dependence of the RHS of (|l^) on T, and e is a 
decreasing function of temperature. This picture is valid 
down to 



at this temperature e saturates becoming 



T K (l + 2|^ln — 



(20) 



(21) 



The quantity Tk gives the low-energy scale of the prob- 
lem and further will be identified with the Kondo tem- 
perature. Note that it has an exponential dependence on 
the order parameter 

The second equation defining the saddle point is ob- 
tained differentiating ( |lC)|) with respect to 4>. There is 
always a trivial solution = 0, the second one is found 
from: 



o 



(22) 



where the summation over discrete ui n was transformed 
to the integration, which is valid in the low temperature 
limit. Assuming e< 1 and eU -C 1, i.e. the system is in 
the vicinity of the phase transition, (E3) becomes: 



1 
J 



dx (l + |0| 2 ln(l + ^))ln(l + 4 r ) 

J * (l+|0|2 1n(l + ^)) 2 +x2[/2 ln 2 (1 + J T) - 



(23) 



The RHS of (23) is bounded from above with the value 
1 achieved at <f) = and U = 0. Therefore we conclude 
that if J is less than the critical value J c given by 



1 



dx ln(l + ^r) 

Jc(U) J 7T l + x 2 U 2 ln 2 (l + \Y 

o x 



(24) 



then (f> — at the saddle point and the impurity is de- 
coupled from quasiparticles. The critical coupling is an 
increasing function of U . In particular 



J C {U)~ l + ^U 
J C (U) » u 



at U < 1, 
at U > 1. 



(25) 



As J increases and becomes larger than J c , the nontriv- 
ial solution of (E23|) becomes relevant, since it defines the 



minimum of the action as a function of \(f>\ 2 . Near the 
critical point J w J c we have: 



j - j c 



i 



r/ 2 



4J2]n2 V 2 In 2 



101 



, J-J C {U) 



Jo(U) 



at {/ < 1, 



at J7 > 1, (26) 



At the transition Tk = and according to ( p0| ) its asymp- 
totics are: 



T K (J, U) « exp (~ a( jg%^ ) at ^ « 1. 
TjcCJ, CO « exp (- 4(J _g ( ( g )ain a ) at U » 1. 



(27) 



These expressions show that close to the phase transi- 
tion Tk is exponentially small and as long as Tk < T <C 
1, the system is in the regime where the effective impu- 
rity energy e is a temperature dependent quantity (see 
1)): 



1 A /-Jc, \ (J-Jc)u TT , 



2(J- J c ([/))lnC7 
J c (t/) + 2( J — J C {U)) In C7 



at 17 > 1. (28) 



The exponential behavior of the Kondo temperature and 
the effective impurity energy has been predicted already 
in p7| . However, we note that this is the case only if the 
particle-hole symmetry is broken. In the particle-hole 
symmetric case e and Tk are identically equal to zero 
for the whole range of the magnetic coupling constant. 
In the regime, where T -C Tk, e saturates at the value 
proportional to Tk (|2l]). 

As we noted in the introduction, our results are quite 
different from those derived in 16 ItJ . We believe there 



are several weak points in those treatments. In particu- 
lar, the action (2.28) in [[Hi) clearly doesn't give the zero 
saddle point for the effective impurity energy e/ at the 
particle- hole symmetric case Qf — N c /2 (we use nota- 
tions of |l6| here). As a result the action in that form 
doesn't predict the Curie law for the susceptibility even 
when the impurity is free. On the other hand this sym- 
metrization is appropriate for finding the order param- 
eter (f> and the critical coupling g c (J c in our notation). 
That is why we find completely different expressions for 
the Kondo temperature as compared to |lql , while the 
critical coupling J c remains the same. In their second 
paper |^7| , the authors considered only the linear den- 
sity of states (see (3.33) therein) and treated the action 
more carefully. However the saddle point equation (5.3) 
in JTtJ is still incorrect. It gives the nonzero value of e/ 
in the particle-hole symmetric case (see e.g. (5.9) with 
x = 1/2). It seems that there was lost a contribution to 
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FIG. 1: Phase diagram in the J—U plane. The solid line (J — 
J c ) corresponds to the Kondo transition at zero temperature. 
At this coupling the order parameter (j> first becomes non-zero. 
The dot and the dash line are the contours of constant Kondo 
temperature Tk{J,U) — 10~ 4 and Tk = 10~ 5 respectively. 
The critical spin-spin coupling J c remains finite at U = 0, 
however in the particle-hole symmetric case there is no real 
transition to the Kondo phase since at any finite temperature 
the crossover value of J (dash and dot lines) diverges as U — > 
0. For r < 1 the phase diagram is similar to this, while for 
r > 1 the zero temperature boundary between the Kondo and 
Curie phases is analogous to the dash line. The solid line in 
this case doesn't represent any real phase transition. 



(5.3) from the residue at |e| > |Ao|, and hence an implicit 
particle-hole asymmetry was introduced. 

Figure 1 illustrates the phase diagram in the J—U 
plane. The solid line corresponds to the phase transition, 
while the dash and dot lines show the contours of the 
constant Tk- At fixed T the latter lines separate high- 
temperature {Tk < T) and the low-temperature {Tk > 
T) regions. Here we would like to emphasize that the 
phase diagram with the dash line b eing a boundary is 
reminiscent to that obtained by NRG [g(J ■ This probably 
indicates that the large N-theory is not very reliable in 
the zero temperature limit. It is also important to note 
that the curve J{U), corresponding to Tk =const is very 
insensitive to the variations of the latter. Thus the dot 
and the dash lines corresponding to Tk — 10~ 4 and Tk — 
10 -5 almost coincide. 



A. Generalization to the nonlinear density of states 

Before we proceed further with D-wave superconduc- 
tors, let us generalize obtained in the large N limit ex- 
pressions to the case of vanishing density of states with 
an arbitrary exponent: 



r+l I 



P(e) = 



|e|<l 
\e\ > 1 



(29) 



I.e. the action is still given by the equation ( |To| ) but the 
unperturbed Green function G[](a; n ) is now as follows: 



T I 1 i+I 3+r _ 1 
.r + l \ ' 2 ' 2 



•(30) 



where T is the hypergeometric function. High and low 
frequency asymptotics of Gq are 



r r+l 7r 
n 2 cos 4 



i^n^ZJ U n < 1 
0J„ > !• 



(31) 



For r < 1/2 the problem was studied in detail in [jTJ], fL5| . 
Both the consequent analysis and numerical results |2C| ] 
show that the effects of the particle-hole asymmetry be- 
come crucial for r > 1/2. Here only this situation will be 
considered. Doing the same steps as before we find the 
critical J, where \4>\ first becomes nonzero, is 



Jc 



r{r + 1) 



J c ~ U max 



2r 
r+l'' 



at U 
at U 



(32) 



In the case of the small potential scattering U, the order 
parameter near the phase transition is: 



J - J c 4r - 2 
J 2 C (r + l) 2 



1> 



1 



'•I 5 



(33) 



where tj) is the digamma function. Note that regardless 
of r, \<j)\ 2 is proportional to J — J c . However, when e is 
considered, it is necessary to distinguish between r < 1 
and r > 1. In the former case one can define the Kondo 
temperature: 



T 



K 



2cosf^ 



4 r_ 1 C0S 2":£(2r) 



_7r r (l-r) 4 r ~! 7r 1 +'-|^| 4 j;((7) 
where £ is the Riemann zeta function, 



(34) 



F r {U)' 



XJ-^- - 1 71-^(1 + r)- 
1 r2 1 r L sin^ 

2r 



for U < 1 
for U > 1. (35) 



(cos^ 

For r < 1 near the phase transition we have: 

T K cx (J -J c )^ (36) 

Above the Kondo temperature e increases with T, while 
for T < Tk it saturates at the value: 



, y+ r {l + r) ( 2cos^ 
2c os^ nr{l + r) K 



(37) 
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It is easy to check that the expressions for r = 1 obtained 
earlier are consistent with these formulas. 

The picture becomes quite different for r > 1. If the 
order parameter is sufficiently small, e can be found from: 



tanh 



/3e 



= 21 



•1, 



I 



W)F r {U). (38) 



Provided the RHS of (|3S]) is less than 1, this equation has 
a solution with e a T, i.e. vanishing as T goes to zero. 
On the other hand if the RHS becomes greater than 1, 
this equation has no solution and at zero temperature at 
this point there is a phase transition to the Kondo phase. 
Thus for r slightly greater than 1 and U < 1 the new 
critical coupling J c is 



J c - J c + 0.35 



1 



U 



(39) 



For (/>lwe have J c ~ J c - The diverging behavior of J c 
at U — > is in fact very similar to that of the crossover 
spin-spin coupling for r = 1 at finite temperature (see 
dash line on fig. 1). However, for r > 1 this divergence 
exists even in the zero temperature limit. 

At this point we get another evidence that the large- 
N theory might be not reliable at small temperatures for 
r < 1. Namely, the NRG analysis shows that r = 1 is not 
a special point in the phase diagram plj . On the other 
hand, the saddle point results predict that the curve sim- 
ilar to the dash (not solid) line in figure 1 represents the 
actual transition to the Kondo phase for r > 1. At the 
same time at finite temperature the mean field phase di- 
agram is qualitatively similar to the NRG result [^5|, . 
And moreover, for r close to 1 the contours Tk(J, U) = T 
are very insensitive to the value of T (see figure 1). 



III. MAGNETIC SUSCEPTIBILITY AND 
QUASIPARTICLE DENSITY OF STATES. 

We start this section from calculating a local magnetic 
susceptibility at zero magnetic field. In the saddle point 
approximation it is given by: 



Xlc 



1 d 2 S eff (h) _ 1 d 2 S eff (e) 



4/3 dh 2 4/3 
Using ( |l2|) we immediately get: 

N 



de 2 



Xlc 



E 



4/3 K-ie-#| 2 G T K)r 



(40) 



(41) 



This expression becomes particularly simple if there is no 
potential scattering and e = (we always assume that the 
temperature is small compared to the cutoff, i.e. 8 3> 1): 



Xlc 



N8 
~16 



l + 2|0| 2 log 



8 



(42) 



'c 

ZJ 

1 

CO 



Q_ 

CD 
O 

o 

ZJ 

CO 

"cO 
Q 
O 




J=2.4 U=0 
J=2.4 11=0.1 
J=2.4 U=0.2 
J=2.4 U=0.25 



20000 40000 60000 80000 100000 

Inverse Temperature (3 



c 

ZJ 

_b 

l 

CO 



Q_ 

CD 
O 
O 
ZJ 

CO 




20000 40000 60000 80000 100000 

Inverse Temperature (3 

FIG. 2: Femperature dependence of the local susceptibility 
(top) and the impurity susceptibility (bottom) in the Kondo 
phase for different potential couplings. For [7 = and U = 0.1 
(see Figure 1) the Kondo temperature is much smaller than 
T = for all plotted values of j3. So both susceptibilities in- 
crease with j3. Note that \ioc is clearly a sublinear function of 
inverse temperature even in the particle-hole symmetric case. 
For U — 0.2 and U = 0.25 the Kondo temperature becomes 
large enough and both susceptibilities have a maximum at 
J3Tk ~ 1. Fhe local susceptibility saturates at j3 — > oo, while 
Ximp — > (see Eq. M). 



It is not surprising that for <f> = in the physical case N = 
2, (^||) gives one half of the usual Curie constant. This 
feature is an artifact of the large N techniques applied 
here. For r < 1 instead of p2) we have 



NB 

Xloc ~-Z- K 8 < fa = 7T 

16 



" 2cos^ 



_7r(r + 1) 



2Np 2r ~ 1 cos 2 2? A r - 1 , s 
^^C(2r) 3» 



(43) 
A)- (44) 



On the other hand for r > 1 the local susceptibility in 
the particle-holc symmetric case is completely described 
by the Curie law. 

The formulas above are valid as long as e <C T. For 
r < 1 this is equivalent to T > Tk- Below the Kondo 
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temperature xio 



Xlc 



Xlc 



saturates at: 
N\^ 



4T K (l-2|0| 2 ln(T K ))3 
N 2(1 - r) cos 2 



K 



for r = 1, (45) 
for r < 1. (46) 



The temperature dependence of Xioc for r = 1 is plot- 
ted on the top graph on Figure 2. This function is non- 
monotonic with the maximum occurring at T w T%. 

Using equation (f4l|), it is not hard to calculate a cor- 
rection to xioc(T = 0) for T <^T K : 



XiociT) - xioc{0) 



9C(3)H 2 T 3 



(47) 



showing that xioc decreases at T — ► for T < Tk- 

The local susceptibility is the response function to the 
magnetic held coupled only to the impurity. We will also 
consider here the impurity susceptibility which addition- 
ally includes the contribution from the nearby electron 
cloud . 

Assuming that g-factors of band and impurity elec- 
trons are the same we see that the only effect of the 
magnetic field h on the action is the shift of the 
imaginary frequency ui n — > u> n + ih for the spin up states 
and u> n — > Lo n — ih for the spin down states. Therefore, 
the impurity susceptibility at zero field is: 



Ximp 



N ^ |0|2^L Gt K) 



43 



E 



If e = 0, expression above reduces to: 

m 2 



Ximp 



N0 
~L6~ 



1 - 



l + 2|^| 2 ln/3 



(48) 



(49) 



Note that Ximp asymptotically approaches Curie law at 
T — > 0, contrary to (42) where Txioc — > 0. This pecu- 
liarity was already predicted in (24] , however in the later 



work |2C)fl , xioc was found to be proportional to the inverse 
temperature without any logarithmic corrections. The 
controversy between the mean-field and NRG results can 
not be resolved considering the next terms in the l/N 
expansion near the saddle point. In fact, it is easy to 
show, that for U — at each order of the perturbation 
series for xioc a multiplier (1 + 2|</>| 2 ln/3) 2 is generated 
in the denominator. Therefore the sublinear dependence 
of xioc{T) is generic for the large iV approach. Prob- 
ably the reason for the discrepancy between this work 
and |2^, |2lJ is that the large N expansion is asymptotic 
in the particle- hole symmetric case and for r < 1, and it 
diverges as temperature goes to zero. 

For r > 1 near J = J c , Ximp = Xioc = N0/16. While 
for r < 1 there is a crossover between 



and 



Ximp 



Ximp 



N0 
~16 



Nr0 

HP 



(50) 



(51) 



T0 O > 1 and T0 O < 1 respectively, O is defined in Q). 
So the impurity is partially screened. We note again that 
these results are valid only for r > 1/2, otherwise there 
is a transition to Kondo phase even in the particle-hole 
symmetric case. 

Below the Kondo temperature Ximp goes to zero in 
agrement with p7| . For example at r = 1: 



T|0| 2 iVln2 
Xlmp ~ r£(l-2M 2 lnT)" 



(52) 



Impurity susceptibility as a function of temperature is 
plotted on the bottom graph on Figure 2. 

Another manifestation of the coupling between quasi- 
particles and impurity is the change of the electron den- 
sity of states, which is given by the imaginary part of the 
effective Green function: 



p(n,r) =QTr 



1 



^G(r,r,fi), 



(53) 



where f2 is the real frequency with an infinitesimal posi- 
tive imaginary part, 



G(r, r', n) = G(r, r', Q) + G(r, r , 0) T K (Q) G(r , r', fi). (54) 
The scattering matrix 7}r(f2), corresponding to the Kondo contribution is equal to: 

Tl<m = -n + er z -\l\l z G{v,,v ^W ^ 

note that a similar form for T was found in p9| j for the case of an Anderson impurity below the Kondo temperature. 
Using (Q), ( |S3"| ) can be simplified further: 

p(Q, r) = 9 G°,(r, fi)T + (0)G° (-r, Q) + G°(r, Q)T_(fi)G£(-r, Q). (56) 
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We adopted the notation: 



G° + (r,Cl) = Tr^^G (r,Cl), G° x (r, CI) = Tr ^ G (r, O), (57) 



T±(fi) i±^(0) + (i±t7Gg(n)J _ 0±e _|^_^_- (58) 



The spatial dependence of DOS coincides with that 
obtained earlier for the case of the pure classical scatter- 
ing pCj ]. We note that at small frequencies (CI <C A), 
there is a maximum at the nearest neighbors of the im- 
purity and at r — > oo the asymptotic behavior of the 
DOS is r~ 2 . Also it is easy to show that if the resonance 
frequency is much smaller than the superconducting gap, 
the spatially integrated DOS remains particle-hole sym- 
metric. It is determined by the momenta in the vicinity 
of the nodes and therefore this result is valid regardless 
the microscopical details of the model such as chemical 
potential, existence of second-neighbor hopping matrix 
elements, etc. 

Clearly the T-matrix for the Kondo scattering channel 
has a pole at the frequency CI ss Tk- If Tk <C A, i.e. 
J ~ J c , the resonant levels corresponding to the poles in 
7+ and 71 become very narrow and they have roughly 
Lorentzian shape. In the opposite limit A <C Tk, the res- 
onance corresponding to 71 disappears and (|56|) results 
in the Fano lineshape in agreement with pq . 

We would like to emphasize that the scattering on a 
strong classical magnet gives two identical peaks at posi- 



tive and negative frequencies and additional weak poten- 
tial coupling results in the weak splitting of those [fl0|| . 
In the Kondo case, even at small U the spatial structure 
of the resonances at positive and negative frequency is 
highly asymmetric. 



IV. NON-MAGNETIC IMPURITY IN A 
D-WAVE SUPERCONDUCTOR 

In this section we focus on the situation, when the 
spin is not located on the impurity site, but rather it is 
a staggered moment distributed nearby. The main mo- 
tivation to this investigation comes from the NMR ex- 
periments H H |7[ |8[ H with non-magnetic Zn or Li 
impurities in a BSSCO. The effective strength of the in- 
duced moments rapidly decreases with the distance ||, 
therefore it is reasonable to assume that the spins are 
sitting only on the first nearest neighbors. 

In this way the new interacting part of the Hamiltonian 
becomes: 



H int = u4(r K j ^(r ) - - ]T(V>}(r s )F T - Fty;(r,))(Ffy T (r.) " (59) 

seo 

where O denotes the subset of the sites closest to the i"o. Contrary to the spin-spin coupling which is strongest on 
the nearest neighbors of the impurity, the potential scattering is dominated by the local term, therefore its form is 
unchanged as compared to (^). Repeating the same steps as before and using the static approximation we obtain: 

S = -Tr\nl^- + er z - £ <j> s G (v s ,r s ,, -pj S , J + M £ (60) 
\ s,s'eo ^ 'I seo 



This action can be also evaluated in the saddle point 
approximation. The important difference with the lo- 
cal spin case is that there are four nontrivial saddle 
points corresponding to S (0io = 0oi = <^ = ^ot)> 
D (</>io = -0oi = 0t,o = ~0o,t)> and p (0m = ~0T,o ! 
0o,i = 0o T = 0) symmetry of the order parameter, the 
latter case being twice degenerate. On the Figure 3 we 
show the dependence of the critical coupling J c on the 



potential scattering for the order parameter of the S and 
D symmetry (P case is always in between) and for the lo- 
cal spin. Clearly, for the whole range of U, D order has a 
lower J c and as a result a smaller action. Also, contrary 
to the local spin case, J c decreases with U. This phe- 
nomenon is readily understood since the potential scat- 
tering increases the quasiparticle density of states on the 
nearest neighbors pQ|, where the magnetic moments are 
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0.2 0.4 0.6 0.8 

Potential Coupling U 



1.0 



FIG. 3: Critical spin-spin coupling as a function of potential 
scattering for the non-local spin case. The D-wave configu- 
ration of the order parameter (j> always has a lower critical 
coupling J c and the smallest action. Contrary to the local 
spin case the Kondo transition occurs at lower value of J 
with increasing U. 



located, and hence enhances the effective spin-spin cou- 
pling. It is interesting to note, that the effective DOS of 
the "bath" electrons for the D-wave saddle point in the 
absence of potential scattering term is proportional to w 3 
and not u> due to Green functions cancellation near the 
nodes |27|] . However this feature is not important for our 
consequent analysis. 

Both local and impurity magnetic susceptibility prac- 
tically don't change as compared to the local spin case, 
since they are sensitive to the value of the impurity spin 
and not to its distribution. The situation is quite dif- 
ferent for the density of states . In particular, for the 
spin scattering part and _D-wave order parameter ( |56| ) 
generalizes to 



where: 



r±(o) 



W(G d (Q)- 



Gi 2 (v a M)Uy°(Q,) 

i-u 2 cf(n) 



(61) 



(62) 



G d +iX (r,n) = - s 2 y )G + , x (r,r s ,tl), 



G d Q (n)= £ (sl-s 2 y )( Sx 2 -s^)G° (r s -r s ,,n), 

s,s'£0 

with Si — 0,±1. The resonance frequency is approxi- 
mately equal to the Kondo temperature as in the local 
spin case. If the potential scattering is not too large 
UG°(T K ) < 1 then the spatial dependence of the DOS 
is determined by the interference from the scattering on 



the four nearest neighbors to the impurity sites. In par- 
ticular, for the D (S) order parameter, the resonance 
corresponding to T_ (T + ) matrix gives the largest DOS 
on the impurity site r = ro, local maximum on the sec- 
ond neighbors and very small signal on the first neigh- 
bors |13| . This behavior is in the excellent agreement 
with the experiments, where the spatial distribution of 
DOS near the Zn impurity was directly observed ||. Also 
at r — > oo the density of states vanishes as r~ 4 , i.e. much 
faster than for the the local spin case. The spatially in- 
tegrated DOS is generically particle-hole asymmetric. In 
particular it can be shown that if UGq(Tk) -C 1 and 
Tk *C A then at the resonance frequency the integrated 
DOS is: 



P+ (Cl) 



A 2 



F 2 A 2 



d 2 k 



4 
A 



(st + Alf W 
A\ d 2 k 



3T+(ft) (63) 



„ , , ,_3r_(-ft) (64) 

I J (e 2 + A 2 ) (2^) 2 



Simple power counting shows that (63) and (|6J) converge 
near the nodes. Therefore the ratio of the spatially in- 
tegrated DOS is nonuniversal and while for particle-hole 
symmetric spectrum p-/p+ ~ 2.5, it decreases with the 
doping crossing 1 at some moment. For comparison we 
will give the results for the local spin: 



p + (Q) 



d 2 k 



£_(-n) 



A? 



A|) W 
d 2 k 



rf+A 2 ,>fi 2 



(4 + A 2 f (2-) : 



QT+(n) (65) 



3fT_(-fi) (66) 



This expressions logarithmically diverge when e and A 
become small, therefore it is necessary to put a low en- 
ergy cutoff of the order of frequency. The integrals mainly 
sit on the nodes and p+ rs p_ if Vt is small. Physically 
this differences are related to the slow (l/r 2 ) dependence 
of the DOS on the distance for a single spin, while for the 
nonlocal scattering the DOS rapidly decreases away from 
the impurity and the main contribution to the spatial in- 
tegrals comes from the short distances or high momenta. 

Let us say a few words about strong potential scatter- 
ing case U 3> 1. As we see from the Figure 3, the critical 
Kondo coupling J c goes to zero, in fact as 1/U up to log- 
arithmic corrections. Therefore for the most values of J 
the system is in the Kondo phase. If J is not too large, 
then the order parameter \<fi\ and the Kondo temperature 
are small enough so that the condition UGq(Tk) <§C 1 ful- 
fills. As a result the spatial structure of the DOS remains 
the same as if there is no potential scattering. This con- 
tradicts the naive idea of the " hard- wall" impurity, which 
prohibits large DOS at r = ro. And the reason, why this 
argument is not correct is that the energy scale Tk is so 
small that the potential scattering is irrelevant. If both 
conditions U > 1 and J > J c are fulfilled then UG°(T K ) 
might be greater than 1. In this case the Kondo scatter- 
ing occurs on the almost frozen density of states created 
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by the potential scatterer and apart from small variations 
the spatial structures of the Kondo and potential reso- 
nances coincide. But we emphasize again that in order 
to have this situation both U and J should be very large. 



V. CONCLUSIONS 

We have studied the screening of a magnetic impurity 
in D-wave superconductors. It is found that there is a 
Kondo transition, which occurs only if the particle-hole 
symmetry is distorted, e.g. by the potential scattering. 
The Kondo temperature is exponentially small near the 
critical coupling in agreement with [ p"7| . However, at fi- 
nite but small temperature the crossover from the regime 
T Tj(, where the impurity spin is not screened to 
T <C Tk occurs at qualitatively different values of the 
magnetic coupling J (sec (fig. 1)), in particular J — > oo 
if the potential coupling U — ► 0. In the particle-hole 
symmetric phase it was found that the impurity suscep- 

const 



imp 



tibility behaves in a Curie-like fashion: T\ 
as T — > zero, while the local susceptibility has multi- 
plicative logarithmic corrections, so that: Txioc — * 0. 
This result, which holds in any order of l/N expansion, 
agrees with [24], however it contradicts later NRG in- 
vestigations [20, [2]]]. It is possible that the logarithmic 
corrections are beyond the accuracy of the large N ex- 
pansion in the particle-hole symmetric case. On the other 
hand, the resulting phase diagram at finite temperature 
is qualitatively similar to that, obtained by NRG [P0[ |21 



and dynamic multichannel P2[ methods. Moreover the 
contours of constant Tk change very slowly (logarithmi- 
cally) with Tk (see fig. 1). Therefore, we expect that 
the conclusions about the Kondo phase are reliable. 

Also we examined the situations with 1/2 < r < 1 and 
r > 1. The former appears to be analogous to the case 
r = 1 with the only difference that exponential (loga- 
rithmical) behavior is substituted by the suitable power 
law. For r > 1 the picture is different in a sense that 
the phase diagram is qualitatively consistent with NRG 
results even in the zero temperature limit. 

Finally we studied the behavior of a non-magnetic im- 
purity in a D-wave superconductor, which creates a stag- 
gered magnetization on the nearest copper neighbors. It 
was found that the saddle point corresponding to the 
true minimum of the action always has a D-symmetry. 
The resulting spatial structure of the quasiparticle den- 
sity of states near the impurity for this saddle point is 
consistent with experimental data. It was also shown 
that unlike the single spin case the critical magnetic cou- 
pling corresponding to the Kondo transition decreases as 
U increases . 
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